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Abstract 

This note provides an interim summary of the current recommendations of the 
PDF4LHC working group for the use of parton distribution functions (PDFs) 
and of PDF uncertainties at the LHC, for cross section and cross section uncer- 
tainty calculations. It also contains a succinct user guide to the computation of 
PDFs, uncertainties and correlations using available PDF sets. 

A companion note (the PDF4LHC Working Group Interim Report) summa- 
rizes predictions for benchmark cross sections at the LHC at NLO using mod- 
em PDFs currently available from 6 PDF fitting groups. 
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1. Introduction 

The LHC experiments are currently producing cross sections from the 7 TeV data, and thus need accurate 
predictions for these cross sections and their uncertainties at NLO and NNLO. Crucial to the predictions 
and their uncertainties are the parton distribution functions (PDFs) obtained from global fits to data from 
deep-inelastic scattering, Drell-Yan and jet data. A number of groups have produced publicly available 
PDFs using different data sets and analysis frameworks. It is one of the charges of the PDF4LHC working 
group to evaluate and understand differences among the PDF sets to be used at the LHC, and to provide a 
protocol for both experimentalists and theorists to use the PDF sets to calculate central cross sections at 
the LHC, as well as to estimate their PDF uncertainty. This note is intended to provide recommendations 
for estimations of cross sections and uncertainties at the LHC. 

2. The PDF4LHC recommendation 

Before the recommendation is presented, it is useful to highlight the differences between two use cases: 
(1) cross sections which have not yet been measured (such as, for example, Higgs production) and (2) 
comparisons to existing cross sections. For the latter, the most useful comparisons should be to the 
predictions using individual PDFs (and their uncertainty bands). Such cross sections have the potential, 
for example, to provide information useful for modification of those PDFs. For the former, in particular 
the cross section predictions in this report, we would like to provide a reliable estimate of the true 
uncertainty, taking into account possible differences between the central values of predictions using 
different PDFs 0. From the results seen it is clear that this uncertainty will be larger than that from any 
single PDF set, but we feel it should not lose all connection to the individual PDF uncertainties (which 
would happen for many processes if the full spread of all PDFs were used), so some compromise is 
proposed. 

The wish for a recommendation follows directly from the HERALHC workshop conclusions HI, 
and has always been one of the main goals of the PDF4LHC group since its creation in 2006, particularly 
as a wish of the LHC experiments. In order for the recommendation to be acceptable by the experiments, 
it has to be pragmatic and not unnecessarily complicated. It is also an advantage to try to keep close to 
the techniques or procedures that are already being used in the experiments up to now. To that end, many 
studies in the past years have been done with CTEQ and MSTW, and recently also NNPDF. But it should 
be clear that at this point no general judgement is made on whether certain PDFs can or cannot be used; 
for any given particular analysis, different expert judgements can lead to different choices, maybe even 
the use of only a single PDF set. Also, the recommendation given below can and will be revised in due 
time when a new level of understanding and development is reached, which is expected to follow from 
the ongoing discussions at the PDF4LHC forum. 

As seen at NLO there is always reasonable agreement between MSTW, CTEQ and NNPDF and 
potentially more deviation with the other sets. In some cases this deviation has at least one potential 
origin, e.g. the tt cross-section at 7 TeV at the LHC probes similar PDFs as probed in the lower-p^- jet 
production at the Tevatron, which has neither been fit nor validated against quantitatively by some groups 
(preliminary results for ABM may be found at [2J). As noted, large deviations in predictions between 
existing NNLO sets are similar to those between the same NLO sets. Discrepancies in MSTW, CTEQ and 
NNPDF do not always have clear origin, or may be a matter of procedure (e.g. gluon parameterisation) 
which is an ongoing debate between groups. Bearing this in mind and having been requested to provide 
a procedure to give a moderately conservative uncertainty, PDF4LHC recommends the following. 

'it may also be more conservative to use this procedure to calculate the uncertainty for the acceptance for a measured cross 
section. 
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2.1 NLO prescription 

At NLO, the recommendation is, for the first case described above, to use predictions from the PDF fits 
from CTEQ, MSTW and NNPDF. These sets all use results from hadron collider experiments, i.e. the 
Tevatron, as well as fixed target experiments and HERA, and they make available specific sets for a vari- 
ety of as{mz) values. The PDFs from these three groups to be used are: CTEQ6.6 H, MSTW2008 B 
and NNPDF2.0 iH. Neither the CTEQ6.6 nor the MSTW2008 PDF sets use the new combined very 
accurate HERA data sets, which are instead used by NNPDF2.0 (updates of the CTEQ (CTIO |6|) and 
MSTW I PDFs will include them) but these PDFs are now most commonly used by the LHC exper- 
iments and are suggested in the recommendation for this reason. The NNPDF2.0 set does not use a 
general-mass variable flavor number scheme (the NNPDF2. 1 PDF set, which does use a general-mass 
variable flavor number scheme is currently being finalized H), but the alternative method which NNPDF 
use for determining PDF uncertainties provide important independent information. 

Other PDF sets, ABKM09 M\M GJR08 HIKni and HERAPDFl.O IM are useful for direct 
comparison to data as suggested for case (2), for cross checks, and for a more extensive and conserva- 
tive evaluation of the PDF uncertainty. For example, HERAPDFl.O allows a study of the theoretical 
uncertainties related to the charm mass treatment. The Os uncertainties (for the PDFs) can be evaluated 
by taking a range of ±0.0012 for 68%c.l. (or ±0.002 for 90% c.l.) from the preferred central value 
for CTEQ and NNPDF. The total PDE-i-a^ uncertainty can then be evaluated by adding the variations in 
PDFs due to as uncertainty in quadrature with the fixed as PDF uncertainty (shown to correctly incor- 
porate correlations in the quadratic error approximation [14]) or, for NNPDF, more efficiently taking a 
gaussian distribution of PDF replicas corresponding to different values of ag lfT5ll20l . For MSTW the 
PDF+as uncertainties should be evaluated using their prescription which better accounts for correlations 
between the PDF and as uncertainties when using the MSTW dynamical tolerance procedure for uncer- 
tainties |16|. Adding the as uncertainty in quadrature for MSTW can be used as a simplification but 
generally gives slightly smaller uncertainties. 

So the prescription for NLO is as follows: 

• For the calculation of uncertainties at the LHC, use the envelope provided by the central values and 
PDF+as errors from the MSTW08, CTEQ6.6 and NNPDF2.0 PDFs, using each group's prescrip- 
tions for combining the two types of errors. We propose this definition of an envelope because the 
deviations between the predictions are as large as their uncertainties. As a central value, use the 
midpoint of this envelope. We recommend that a 68%c.l. uncertainty envelope be calculated and 
the as variation suggested is consistent with this. Note that the CTEQ6.6 set has uncertainties and 
as variations provided only at 90%c.l. and thus their uncertainties should be reduced by a factor 
of 1.645 for 68%c.l.. Within the quadratic approximation, this procedure is completely correct. 

2.2 NNLO prescription 

At NNLO, base the calculation of PDF uncertainties on the only NNLO set which currently includes a 
wide variety of hadron collider data sets, i.e. MSTW2008B There seems to be no reason to believe that 
the spread in predictions of the global fits, i.e. MSTW, CTEQ and NNPDF, will diminish significantly 
at NNLO compared to NLO, where this spread was somewhat bigger than the uncertainty from each 
single group. Hence, at NNLO the uncertainty obtained from MSTW alone should be expanded to 
some degree. It seems most appropriate to do this by multiplying the MSTW uncertainty at NNLO by 
the factor obtained by dividing the full uncertainty obtained from the envelope of MSTW, CTEQ and 

^The MSTW presentation at the DIS 2010 workshop (7} can be consuUed to assess the effects of these data. 
^The NNPDF presentation at the DIS 2010 workshop | 8| can be consulted to assess the effects of these corrections. 
''Although inclusive jet data from the Tevatron are included in the MSTW2008 (and other) NNLO fits, we note that to date 
the inclusive jet cross section has only been cslculated to NLO. 
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NNPDF results at NLO by the MSTW uncertainty at NLO. In all cases the as uncertainty should be 
included. We note that in most cases so far examined for the LHC running at 7TeV centre of mass 
energy this factor of the envelope divided by the MSTW uncertainty is quite close to 2, and this factor 
can be used as a short-hand prescription. 

Since there are NNLO PDFs obtained from fits by the ABM,GJrEI and HERAPDF groups, these 
should ideally be compared with the above procedure. 

So the prescription at NNLO is: 

• As a central value, use the MSTW08 prediction. As an uncertainty, take the same percentage 
uncertainty on this NNLO prediction as found using the NLO uncertainty prescription given above. 



^See Ref. ||17il for a comparison of a number of benchmark cross sections at NNLO. 
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3. The PDF4LHC prescription for Higgs production via gluon fusion 

The total hadronic cross-section for the inclusive production of a Higgs boson via gluon fusion can be 
written as 

a{hih2 H + X) = [ dxidx2 fa,hA^^^I^F)fb,h2i^'^^f^F) [ dz 6 ( z - aab{z) (1) 

Jo \ X1X2J 

The accurate evaluation of the partonic cross section dab is the most important step to obtain an 
accurate prediction of the central value of the total hadronic cross section. In fact the latter receives large 
contributions from NLO- and NNLO-QCD corrections, but also from the soft gluon resummation, from 
leading NNNLO-QCD terms and from NLO-EW corrections. Finite top mass corrections at NNLO- 
QCD have been studied and turn out to be small. For the sake of simplicity but still in full generality, 
the evaluation of the size of the combined PDF-i-as uncertainty can be obtained by considering only exact 
NLO-QCD and NNLO-QCD corrections in the infinite top mass limit. The results in the next Sections 
have been computed using the code described in [18], improved with the NNLO-QCD conections [19J. 

In this study |^ the following three global sets of PDF have been considered: CTEQ6.6 IH, 
MSTW2008 Elllll, NNPDF2.0 |5|. The combined PDF-t-a^ uncertainties for each of the three global 
sets is computed as discussed in the PDF4LHC Working Group Interim Report. The various recipes are 
also compared in detail in Ref. |15j. 

In order to obtain a meaningful comparison of the predictions computed with different PDF sets, 
it is crucial to adopt the same uncertainty range for the value. For the purposes of illustration of the 
PDF4LHC recommendation, we will assume the same range as the one used for the PDF4LHC Working 
Group Interim Report namely, 

<5(90)a^ = 0.002 90 %c.l., 5^^^^as = 0.0012 = O.OO2/C90 68 %c.l. (2) 

where C90 = 1.64485 is the number of standard deviations corresponding to a 90% c.l. 
At NLO-QCD, the PDF4LHC recipe can be summarized in the following steps. 

1. Compute the Higgs cross section, using CTEQ6.6, MSTW2008nlo68cl, NNPDF2.0. 

2. For each set, use the preferred Osimz) value (respectively 0.118, 0.1207, 0.119). 

3. Compute the PDF-i-as uncertainty band, according to the rules of each collaboration described in 
the PDF4LHC Working Group Interim Report 

4. Take the envelope of the three uncertainty bands. 

5. Compute the mid-point of the resulting band and call uncertainty the distance of the edge of the 
envelope from it. 

In Fig. [T] we show at different collider energies (Tevatron, LHC 7 and 14 TeV) the size of the 
combined PDF-i-a^ uncertainty bands obtained with CTEQ6.6, MSTW2008nlo68cl, NNPDF2.0, all nor- 
malized to the central value by MSTW2008nlo68cl. One remarks that for different Higgs mass values 
the predictions show partial agreement of different pairs of the three collaborations, in such a way that 
only an envelope (the black line) of the three bands provides a conservative estimate of the uncertainty 
in the whole mass spectrum. 

Now we turn to discuss NNLO-QCD. There is, at present, only one global set of PDF extracted 
with NNLO-QCD accuracy: MSTW2008nnlo. In this case it is not possible to prepare an envelope 
like it is done at NLO-QCD. The PDF4LHC recipe in this case extrapolates the information available at 
NLO-QCD to NNLO and can be summarized in the following steps. 

1. From the NLO-QCD exercise determine the percentage uncertainty relative to the mid point of the 

envelope. 

*We would like to thank Alessandro Vicini for carrying out this study. 
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2. Determine the ratio of this percentage uncertainty to the one estimated with MSTW2008nlo68cI. 

3. Compute the central value and the PDF+a<jUncertainty band with MSTW2008nnIo68cI. 

4. Rescale the MSTW2008nnIo percentage uncertainty by the above ratio. 

One remarks from Fig.|2]that the PDF+as uncertainty bands obtained with MSTW2008 at NLO- and at 
NNLO-QCD are very similar. The small differences are taken into account by rescaling the MSTW2008nnlo68cl 
uncertainty band with the ratio at NLO-QCD of the percentage width of the envelope with respect to its 
mid point, to the percentage uncertainty of MSTW2008nIo68cl. As it can be observed in Fig. [2] (lower 
panel), the rescaling factor is close to 2, but it has a non trivial dependence on the Higgs boson mass, on 
the collider type and on the collider energy. 
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Fig. 1: Combined PDF+QsUncertainty band relative to the total Higgs production cross section via gluon fusion, at NLO-QCD, 
evaluated according to the PDF4LHC recipe. The bands are normalized with the central value of MSTW2008nlo. 
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Fig. 2: In the upper panels, combined PDF+asUncertainty band relative to the total Higgs production cross section via gluon 
fusion, at NLO-QCD and at NNLO-QCD, obtained with MSTW2008nlo68cl and with MSTW2008nnlo68cl. In the lower 
panel, rescaling factor obtained from the ratio of the percentage width of the NLO-QCD envelope with respect to its mid point 
to the percentage uncertainty of the MSTW2008nlo68cl band. The rescaling factor has to be applied to the MSTW2008nnlo 
PDF+Qsband, to obtain an estimate of the NNLO-QCD envelope. 
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4. Summary 

In this note, we have provided a method to calculate reasonable estimates for PDF uncertainties for cross 
section predictions at the LHC, using PDFs from CTEQ, MSTW and NNPDF. The method is intended to 
be both pragmatic and practical for the calculation of PDF uncertainties; it is not intended to discourage 
comparisons to all PDFs relevant for LHC calculations. 

The recommendation is expected to evolve when new experimental sets and new PDF determina- 
tions warrant. 
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